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ABSTRACT 

\ The observation of g-mode candidates by the SoHO mission opens the possibility of probing 

^ ■ the internal structure of the solar radiative zone (RZ) and the solar core more directly than 

OO ' possible via the use of the p-mode helioseismology data. We study the effect of rotation and 

(N . RZ magnetic fields on g-mode frequencies. Using a self-consistent static MHD magnetic field 

model we show that a 1 % g-mode frequency shift with respect to the Solar Seismic Model 
(SSeM) prediction, currently hinted in the GOLF data, can be obtained for magnetic fields as 
\^ \ low as 300 kG, for current measured modes of radial order n = — 20. On the other hand, we 

■ also argue that a similar shift for the case of the low order g-mode candidate (I — 2, n = —3) 
"^H I frequencies can not result from rotation effects nor from central magnetic fields, unless these 

p exceed 8 MG. 

O ' Key words: MHD - Sun: helioseismology - Sun: interior - Sun: magnetic fields. 
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1 INTRODUCTION 



The observation of g-mode candidates by th e Global Oscillation at L ow Frequencies (GOLF) instrument aboard the ESA/NASA Solar and 
Heliospheric Observatory (SoHO) mission tTurck-Chieze et al .|2004) opens new prospects for solar physics. Indeed, it brings the possibility 



of probing the internal structure of the solar radiative zone (RZ) and the solar core more efficiently and more directly than possible via the 
use of the p-mode helioseismology data. 



The shift of the g-mode candidate frequ encies with respect to the Solar Seismic Model (SSeM) prediction jTurck-Chieze et al.ll2001 



Couvidat. Turck-Chieze & KosovichevI 120030 hinted in the experiment is of the order Sujnim /i^ni ^ !%• Other standard models lead to 
the same kind of discrepancy but it is clear that none of these models take into account the solar internal dynamical effects which might 
contribute to this shift. So one may attempt to explain this either as resulting from a strong central magnetic field or by the rotation of the 
RZ, or both. Here we study the effect of both RZ magnetic fields and rotation on g-mode frequencies in order to check previous estimates of 

their values 0. 

One approach is provided by linearized one-dimensional (1-D) magneto-hydrodynamics (MHD) jBurgess et al .l2004ah . This enables us 



to determine analytically the MHD g-mode spectra beyond the JWKB approximation and opens the possibility of Alfve n or slow resonances 
that co uld produce sizeable matter density perturbations in the RZ, and potentially affect solar neutrino propagation i Burgess et alj 200 
2004bh . 



Magnetic s hifts of g-mode frequencies were calculated using exact eigenfunctions in 1-D MHD, within the perturbative regime (small 
magnetic field) ( Rashb a. Semikoz & Valle 2W(3). However such approach can not describe splittings in g-mode frequencies, intrinsically 



associated with spherical geometry. Moreover, all shifts have the same sign, in contrast with first indications from experiment. 
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^ The resulting magnetic and rotation frequency splittings behave as even (in the Jeffreys-Wentzel-Kramers-Brillouin (JWKB) approximation used below) 
and odd with respect to the azimuthal number m, as seen in eqs. Uli and U2\i . 
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In this paper we generalize the MHD picture to t he three-dimensional (3-D) case using standard perturbative appr oach for the calculation 
of magnetic field corrections to the g-mode spectra ( Unno et al. 19891 : Hasan. Zahn & Christensen-Dalsgaard 2005 ). The main ingredients 
of our calculation are: a) a 3-D model of the background magnetic field in RZ and b) the radial profiles of the eigenfunctions for the hori- 
zontal and radial displacements £,h,rir). F or (a) we will use for t he first time the static 3-D background magnetic field solution introduced 
by Kutvitskii & Solov'ev (KS, for short) i Kutvitskii &_Solov^ev[ lj94). Fo r (b) we will adop t, as a first approximation, the eigenfunction 
calc ulated for the SSeM neglectin g magnetic fields jMathur. Turck-Chieze & Couvidal 200a) using the Aarhus adiabatic oscillation pack- 
age ( IChristensen-Dalsgaaral2003h . 

We organize our presentation as follows. In order to calculate the 3-D magnetic shift we first derive in Se ction 2 the perturbative magn etic 
field B' using the explicit 3-D form of the background magnetic field B obtained in a static MHD model i Kutvitskii & Solov'evlll994 ). In 
Section 3 we derive a simple formula for the magnetic shift in the JWKB approximation (hig h order g-modes, \n\ ^ 1) and find agreement 
with calculations performed for the case of slowly pulsating B-s tars ([Ha san. Zahn & Christ ensen-Dalsgaard 2005 ). The relevant integrals 
were calcul ated numerically using KS background magnetic field {Kutvitskii & Solov'evlll994l) and the eigenfunctions ^^.h were taken from 
the SSeM ( jCouvidat. Turck-Chieze & Kosovichevll2003t) model 2). In Section 4 we calculate the rotational splittings of g-modes in order to 
compare with those induced by the magnetic field. Finally, in the Discussion we summarize our results. We conclude that future more precise 
g-mode observations will open prospects for narrowing down the i nferred magnetic field strength ranges considerably, when compared with 



values currently discussed in the literature, 1 G < i? < 30 MG i Couvidat. Turck-Chieze & Kosovichevll2003l : iTurck-Chieze et alJl2001 



MosJ2003l : iRuzmaikin & Lindsevll2002l : iRashba. Semikoz & ValleH2006h 



2 MAGNETIC FREQUENCY SPLITTING OF G MODES 

We u se the following standard formula for 3-D MHD perturbative correction (see, e. g. i Urmo et al J 19891 : Hasan. Zahn & Christensen-Dalsgaard 
2005h ): 

Slu _ 1 /-[(V X B') X B + (V X B) X B'] dV 



where the eigenvector ^ is given by 
I smO 



(1) 



(2) 



The perturbative magnetic field B' obtained from Faraday's equation in ideal MHD, B = V x (^x B), takes the form: 
B' = [(BV)C, - {£y)Br - Bru] e, + 



+ 



(BV)e. - (fV)Be + (Bo^^-^^B^) _ Ben 



eg + 



where the compressibility u = div ^ is given by 



(3) 



" = + "^"^ ^ — "I s« ^ — ■ 

or r r ad r rsmo ocp 

Ins tead of using an ad hoc an satz, we use the static 3-D configuration for the background magnetic field B in the quiet Sun obtained 
in I Kutvitskii & Solov' evil 1994 ). It expresses the equilibrium between the pressure force, the Lorentz force and the gravitational force. 



Vp - + /9V$ = 0. 

c 

This axisymmetric field 

Ji{r,6) = Bcore {cos9b^{x),sm9b'fj{x),sm9b^{x)^ 



(4) 
(5) 



is specified as a family of solutions that depend on the roots of the spherical Bessel functions labeled by Zk - These come out from /5/2 = 
^fzJ^I-iiz) imposing the boundary condition that B vanishes on the solar surface, x — t/Rq — 1, B(_R n) — 0. The amplitude Bco re 
characterizes the central magnetic field strength. The radial and transversal components of magnetic field are i Kutvitskii & Solov' evlll994 ): 



br{x) = 



blix) 



bl{x) 



(1 - Zk/ sin Zk) 
1 

2(1 — Zk/s'mzk) 
1 

2(1 - Zk/sinzk) 



x^Zk sm Zk 
3 



sin{zkx) 



2 + 



ZkX 

sm(zkx) 



cos{zkx) 



ZkX - 



x-^Zksmzk \ ZkX 
3 / sm(2fcxj 



cos(2;fea::) — ZkXsm{zkx) 



X sm Zk 



ZkX 



COSiZkX^ 



(6) 
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Sun ^ ' *^Sun 

Figure 1. Perpendicular (left), = -^(feg)^ + (fe^)^, and radial (right panel), 6*, components of modes A; = 1, 3, 10 of the self-consistent B-field Eq. (S)- 



Notice that the behavior of B at the solar center (r = 0): 

Br{0,e) = BcoreCOsO, 

Be{0,e) = -Bear, sin 6, 

B40,e) = B,oresm9^-^ ^ 0, (7) 
2 it© 

is completely regular, determined only by the single parameter Bcore- 

In the left panel of Fig. [T] we display the perpendicular component b']_{r) — + (b^)^ while in the right panel we show the 

radial fe5;(x)-component given by Eq. ^ for the different modes k. Note that such low order mo d es of the large-scale magnetic field Eq. l|5j, 
A; = 1, 2, ... 10, survive against magnetic diffusion over the solar age, ~ 4.6 Gyr [Miranda et alJlOOll) . 

Substituting the background magnetic field in Eq. lO into Eq. l|3) one gets all components of the perturbative field B {r,9,(f)), as 
follows: 



Bj, 



Bcore / sin 9deYr 
Rq \ X 



imh^{x) ^^^^ I dbr{x) ^ 2br(x) 



dx 



l{l + l)coseYr ^ , , , 

+ —br{x)^h{x) 



B, 



Rq 

-sin err 



COS OdgYj" 



bflCr I d{bg{x)Ux)) 

X dx 



dx 
+ bg sin ( 



irndgYT, , 



^ (d^vr + i{i + i)Yr) 

X ^ ' 



7-,' Bcor& ( . , fi-^rn 

Bs = — — I i-m cot oYi 



Rq 

-sinerr" 



br^-^{2bg+b,.) 
OX X 



d{^r{x)b4,{x)) ^ irb^ 

dx X 



+ 



sin 6h^^h 



+ im — —dgYi — 

X 

1 d'^Y 



(8) 



sin^ 9 dcfP 

Using displacement eigenfunctions (x), ^J;"'' {x) for the g-modes of the order n and degree 1 1 Mathur. Turck-Chieze & Couvidalll2006h 
we are able to calculate the perturbative magnetic field B (r, 9, 0) given by Eq. ([Sjl and then the integrals in the initial Eq. (Til. This way we 
obtain the magnetic shift of eigenfrequencies. 



3 JWKB APPROXIMATION 



In order to simplify Eq. (|T} we adopt the JWKB well-known approximation llChristensen-Dalsgaardll2003h (see eqs. (7.129)-(7.131) there). 
This holds for large orders n 3> 1 and correspondingly u <Si N (N is the Brunt- Vaisala frequency). The resulting g-mode eigenfunctions 
Ci"''' ^r"'' of oscillations within RZ are given by 
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4 



(9) 
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Figure 2. Horizontal and radial eigenfunctions as a function of normalized radius, r / Rq, are shown for modes gjg (left) and g^^ (right) panel. 
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(10) 



where L = •y/I(T+l) and the JWKB connection d^r/dr ~ L'^£,h/r was used in obtaining ^h. 



Obviously, the ratio between the amplitudes of the horizontal and vertical displacements obeys (see Eq. (7. 132) in dChristensen-Dalsgaard 



20031) ) 



1/2 



showing explicitly that for low degrees 1 — 1,2 the g-mode oscillation is predominantly in the horizontal direction, 2> ^r- Moreover, 
differentiating the horizontal component once more, d^h/dr, one finds that d^h/dr ^ ih/r, since the additional large factor [{N/luY — 
]^ji/2 ^ ^ appears. 

We see that the large derivatives dS^h jdx enter only in the ang ular components B(, and given by Eg. l[p , which with the angular 
factors (cos QdeY^ and (cot QY~) respectively. As a result, as in Ref. i Hasan. Zahn & Christensen-Dalsgaard 20051) . the dominant term in the 
integrand of the numerator in Eq. (TJ is the first one 



-[(V X B ) X B] -C* 



B 



Be, 



l^Afyxhrih) 

X 



da; 



^dY 



+ loot 9Y\ 



The second term in the numerator in Eq. lUll is negligible since it is linear in d^h/ds (in contrast the first term is quadratic). Finally, the 
leading contribution to the denominator is since 2> for ^ 2> 1. 

Under the above approximations we finally obtain the magnetic shift of g-mode frequencies in the Sun, as 



5lO C,„(Bcore)' / h d(xbr^h) / dx\ X^ Ax 



Sttlj^PcR; 



/ -Bcore \ 

V MG / 



(11) 



Notice that, in contrast to the 1-D case, here there is a dependence on the angular degree / and azimuthal number m. This induces a 
splittin g (not only a shift) of the g -mode frequencies. Note also that this expression is in full agreement with the magnetic shift derived 
in ( .Hasan. Zahn & Christensen-Dalsgaard 2005) for the g-mode spectra in slowly pulsating B-stars. Moreover, the coefficient Cim corre- 
sponding to the angular integration coincides exactly with those of Ref. i Hasan. Zahn & Christensen-Dalsgaarj|2005l). Notice however, tha t 
this coincidence was not expected a priori since here we adopt a self-consist ent background magnetic field ( iKutvitskii & Solov'evlll994b 
instead of an arbitrary ansatz for b{x), as in Ref. (.Hasan. Zahn & Christensen-Dalsgaard 2005.) . 



In order to compute the magnetic shift we perform the numerical integration in Eq. dllb using the SSeM eigenfunctions ( IMathur. Turck-Chieze & 



9io 



shown in Fig.|2] Our results are given in Table 1. These results for the magnetic shifts of the g-mode frequencies in the Sun for 
and 320- One finds that magnetic shifts of the magnitude hinted by the data (order 1%) can be obtained for RZ magnetic field strength 300 
kG for the mode 320 ^rid 2500 kG for the mode gjn- We also give the correspon ding values of the Sc coefficient in Eq. Jllb . Note that the 
case of young B-star considered in dHasan. Zahn & Christensen- Dalsgaardliooi) one obtains 110 kG and 11 00 kG respectively. 



Note that for low order g-mode candidates, e.g. gf, (n = —3, 1 = 2, i/^ = 222.02 /iHz) the JWKB approximation is not applicable and, 
moreover, the RZ magnetic field strengths required to provide the magnetic shift ~ 1% would be huge, hence outside our assumptions. 



2 Note we use the exact eigenfunctions, not the approximations f ^ j, given by Eqs. j9), llOt . 
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Table 1. Solar RZ magnetic field which would produce a frequency of 1 % for g-modes g}^^ and g . The magnetic field mode A; = 1 is assumed. See text. 



Mode 


u (pRz) 


P (hour) 


Sc (MG-2) 


Score (MG) 




33.19 


8.37 


1.25 X 10-1 


0.28 


9lQ 


63.22 


4.39 


1.63 X 10-3 


2.48 



Table 2. Here we assumed f! = 430 nHz. 









JWKB, 3D 


3D 


Mode 


V (/iHz) 


P (hour) 












xm 10-3 


xm 10-3 


9\ 


263.07 


1.05 


0.82 


0.95 


9l 


222.46 


1.25 


1.61 


1.54 


9\ 


128.34 


2.17 


1.68 


1.64 


9\o 


63.22 


4.39 


3.40 


3.63 


9w 


103.29 


2.69 


3.47 


3.57 


920 


33.19 


8.37 


6.48 


6.60 


920 


55.98 


4.96 


6.40 


6.44 



4 ROTATION-INDUCED SPLITTING OF G-MODES 

Here we give the standard estimate of the rotational splitting of the g-mode frequencies obtained in the absence of magnetic field i Christensen-Dalsgaarj 



2003h . The exact result reads. 



dUJnlm = m-^ -Tj . (12) 

where the last two terms in the numerator come from the Coriolis force contribution. 

For high order g-modes, n 2> 1, for which JWKB approximation is valid, £,h 3> ^r, and assuming rigid rotation of the RZ, Q = 
constant, one obtains the simple 3D formula: 

SiOnlm = m/J„iQ, (13) 

where the last term in the coefficient P„i ~ 1 — [1(1 + 1)]^^ is given by the Coriolis force contribution. 

For the lowest degree I = 1 this gives P„i = 1/2 while for a high degree I 2> 1 the Coriolis term is negligible and P„i = 1- Choosing 
f2/27r = 430 nHz (T = 27 d) we see from the Table|2]that the JWKB approximation given by Eq. l ll3t (fourth column) is in good agreement 
with the exact calculations using Eq. il2l and given in the last column. 



5 DISCUSSION 



Recent observations of g-mode candidates in the GOLF experiment JTurck-Chieze et al J 1 20040 indicate shifts in the g-mode frequencies 



corresponding to low order modes that may be as large as ~ 1%, e. g. gl, g|. For such modes the expected rotation splittings seem rather 
small. We have found that (5w(rotation) /oj ~ m x 0.001 — 0.0015 for a flat rotation in the radiative zone. As result rotation cannot provide 
the frequency shift hinted by GOLF observations for these modes. Nevertheless as the rotation in the core can be slightly increased by a factor 
greater than 2, it remains important to properly determine the centroid of the observed modes. Moreover, it is important to go beyond the 
standard model picture and introduce the effect of transport and mixing in the solar radiative zone which may also impact on the prediction 
of the high frequency gravity modes. 

Here we mainly study that such shifts may result from RZ magnetic fields. In order to give an estimate of the corresponding values 
of RZ magnetic fields we have adopted the standard 3-D MHD JWKB perturbative approach, aware that its validity is limited to relatively 
large modes. We have found that the magnetic shift of g-mode frequencies for high order modes Slu/u ~ B'^^^^juP' can be of order 
1% for realistic central magnetic fields values obeying all known upper bounds discussed in the literature. For example, we have obtained 
magnetic fields Bcore — 300 kG and Bcore — 2500 kG for 320 ^rid glo, respectively. Such values are in qualitative agreement with 
recent c alculated magnetic field values 110 kG and HOP kG providing magnetic shift 1% for the same modes in the case of young 4 Mq 
B- stars jHasan. Zahn & Christensen-Dalsgaar 31200 5*). Future observations of such high order g-modes n ^ 1 in the Sun could shed light on 



the validity of our assumption. In fact we have already some signature of these modes obtained from their asymptotic properti es. However, 



due t o their very low amplitudes, we have presently mainly analyzed the sum of dipole gravity modes from n = —4 to —26 jGarcia et al 
2006h . A proper decomposition of the observed waves will place more constraints on the shift of these modes. 
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In contrast, to consider the case of low order modes hinted by the recent GOLF data requires further assumptions. In order to provide 
estimates for this case we consider two approximations. First we recall the general dependence Slo/uj ^ B^ore/'^^- It implies that, in order 
to provide the same magnetic shift ~ 1 % the magnetic field for modes q|, should be much str onger than in the case of JWKB frequencies, 
n 2> 1. E.g. for the g-mode candidate (?| with v — 222.46 /iHz ( ITurck-Chieze et al J 120041) from the frequency ratio 1^(53 )/i^(5io) ~ 
222.46/63.22 — 3.52, one finds that the magnetic field which provid es a shift ^1% for would be ~ 8.8 MG. Such large RZ magnetic 
field estimate also agrees with the 1-D MHD approach given in Ref. / Rashba. Semikoz & Valle 2006h (although such 1-D picture can not 
provide angular splittings, it has th e merit of not relying on the JW KB approximation). Indeed, such trend for low order modes, n — 
—3, —2, — i, IS seen m Fig. 4 of Ref. jRashba. Semikoz & Vallel2006h taking into account that the transversal wave number parameter k = 2 
of the 1-D approach is analogous to the degree I in the 3-D approach and 2.8 times more magnetic field value for the maximum Brunt- 
Va isala frequency approximated as A'^ — 2.8 x 10^'^ rad/sec = constant instead of ~ 10^'^ rad/sec applied in numerical calculations 
in jRashba. Semikoz & Vaii3l2006h . 

The possible existence of such strong magnetic fields in the RZ is somewhat disturbing. It could be that neither rotation, nor magnetic 
fields are responsible for the hinted frequency shift with respect to SSeM prediction. For example, the mode might mix with other gravity 
modes penetrating RZ with frequencies just below Brunt-Vaisala frequency. 

In this paper we have considered the possibility of a 1% magnetic shift of eigenfrequenci es, independently of the frequen cy considered. 
The sensitivity to the detailed physics depends on the order of the mode as shown in cMathur. Turck-Chieze & Couvidalll2006h and is greater 
at high frequency (low order) than at large order, demonstrating the great interest of these modes. 

The effect of the magnetic field on the frequency shift is expected to depend on its shape. In a forthcoming paper we plan to explore 
the sensitivity of magnetic shifts of g-modes in the Sun to the structure of magnetic fields given by the self-consistent MHD model. Indeed, 
refined g-mode data may enable future tomography studies of the structure of RZ magnetic fields. 

In short, the search for fundamental solar characteristics deep within radiative zone constitutes an important challenge for helioseis- 
mology methods. Here we have given a conservative limit on the magnitude of the magnetic field obtained from its possible effect on 
the g-mode frequencies. Alternatively, the propagation of neutrinos may probe sh ort-wave MHD density perturbations in the solar RZ 
(\mhd ~ 100 - 200 km) as already discussed in Refs. jBurgess et al.ll2003ll2004j|bl) . 
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